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ABSTRACT 

We present SUSY_FLAVOR — a Fortran 77 program that calculates important leptonic and 
semi-leptonic low-energy observables in the general i?-parity conserving MSSM. For a set 
of input MSSM parameters, the code gives predictions for the K^K^, DD, B^Bd and SgBg 
mixing parameters; B — >■ ^^7, Bs^d l~^l~, K^ — )■ tt^i^z^ and K~^ — j- n'^uu decay branching 
ratios; and the electric dipole moments of the leptons and the neutron. All these quantities 
are calculated at one-loop level (with some higher-order QCD corrections included) in the 
exact sfermion mass eigenbasis, without resorting to mass insertion approximations. The 
program can be obtained from jhttp : / / www . f uw . edu . pi/ susy_f lavor 



1 Introduction 



Flavor changing neutral currents (FCNCs) in the Minimal Supersymmetric Standard Model 
(MSSM) [T] originate from the fact that one cannot, in general, simultaneously diagonal- 
ize the mass matrices of fermions and their supersymmetric partners. The misalignment 
between these mass matrices leads to FCNCs at tree level. Moreover a large number of 
the MSSM parameters which can take complex values is a potential source of CP violation. 
Thus supersymmetric contributions to amplitudes of processes violating flavor and to quan- 
tities measuring CP violation, like Electric Dipole Moments (EDMs), can exceed by orders 
of magnitudes the ones of the SM particles. Such large effects are ruled out by experimen- 
tal measurements which generally agree with the SM predictions and thus provide strong 
bounds on the amount of flavor and CP violation in the MSSM. For instance, measure- 
ments of the kaon system properties prohibit FCNC couplings between the first and second 
generation of down-type squarks larger than 10"'^. These strong limits are often called the 
"SUSY flavor problem" . However, there are also areas where current experiments still leave 
room for large SUSY contributions. For example, constraints from i?-meson experiments 
allow an 0{1) mixing between the second and third generation down-type squarks. Such a 
large mixing could produce FCNC effects that could be observed in the future at _B-factories 
and/or hadron colliders, like the Tevatron and LHC [2]. 

Even assuming the so-called Minimal Flavor Violation (MFV) scenario in which all FCNC 
effects originate from the superpotential Yukawa couplings, the flavor conserving soft SUSY 
breaking parameters can still contain complex phases that cannot be absorbed by a redef- 
inition of fields and can, for example, give large contributions to the electron and neutron 
EDMs. 

As the accuracy of rare decay experiments improves, it is important to have a universal 
computational tool which would help compare new data with the predictions of the MSSM. 
Constructing such a tool is a non-trivial task because finding SUSY contribution to each 
rare decay requires tedious calculations, especially when one wishes to have fully general 
formulae that do not rely on the restrictive assumptions of the MFV scenario. Numer- 
ous analyses have been published in the literature, but because of the complexity of the 
problem, they mostly take into account one, or at most few, rare decays simultaneously. 
Furthermore, most analyses done for general flavor violation in the MSSM soft terms use 
the mass insertion approximation (MIA) (see e.g. [3]) which significantly simplifies calcula- 
tions but does not produce correct results when flavor violation in the superpartner sector 
becomes strong. 

In a series of papers since 1997 PHTB], many supersymmetric FCNC and CP- violating 
observables were analyzed with loop-level accuracy within the setup of the fully general im- 
parity conserving MSSM without resorting to any MIA-type expansions. A FORTRAN 77 
computer programs based on the common set of Feynman rules of Ref. [13] were developed 
for each process. Because these programs use the same conventions, input parameters, and 
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Observable 


Experiment 


AF = 


\de\{ecm) 
\d^\{ecm) 
\dr\{ecm) 
\dn\iecm) 


< 1.6 X 10-27 p] 

< 2.8 X 10-19 [I7j 

< 1.1 X 10-17 |l8j 

< 2.9 X 10-26 [19j 


AF = 1 


Bt{Kl ^ TT^UU) 

Br{K+ 7T+uu) 
BT{Bd ^ ee) 
Br(5d ^ fifi) 
Br(fid ^ rr) 
Br(fi, ^ ee) 
Br(fi, ^ ^;u) 
Br(5, ^ rr) 
Br(E -> X,7) 


< 6.7 X 10-^ [20J 
17.3t\li X 10"^^ [21J 

< 1.13 X 10-7 [22] 

< 1.8 X 10-*^ [23] 

< 4.1 X 10-3 [24] 

< 7.0 X 10"^ [25] 

< 5.8 X 10-*^ |23] 

(3.52 ±0.25) X 10-^ [26] 


AF = 2 


AMk 
AMd 

AMb, 


(2.229 ±0.010) X 10-3 [18J 
(5.292 ± 0.009) X 10-^ ps-^ [18j 

(2.37ir7?) X 10-2 ps-i m 

(0.507 ±0.005) ps-i [26] 
(17.77 ±0.12) ps-i 127J 



Table 1: List of observables calculated by SUSY_FLAVOR and their currently measured values 
or 95% CL bounds (except for Br(i?rf — > ee) and Br(i?d tt) for which the 90% CL 
bounds are given.). 

internal data structures, they can naturally interface with one another. Combining these 
works, we present in this article SUSY_FLAVOR - a publicly available computer code that 
simultaneously calculates the set of important AF = 0, 1, 2 FCNC and CPV observables in 
the framework of the general MSSM. The current version (1.0) of the program takes a set 
of MSSM parameters and calculates the processed listed in Table [U 

Several programs allowing to analyze various aspects of the MSSM flavor phenomenology 
have been published. The most relevant to SUSY_FLAVOR are: CPsuperH [28], SusyBSG [29] 
and Super I so [30]. SusyBSG is dedicated to high-precision predictions for B ^ s'y while 
CPsuperH and Super I so calculate processes similar to the ones computed by SUSY_FLAVOR. 

^In the current version SUSY_FLAVOR calculates also the full one-loop corrections to lepton flavor violating 
i3-meson decays such as _B — > fir. However, it is known that contributions to the amplitudes of these 
processes are greatly enhanced at large tan /3 by formally two loop double penguin diagrams |15j which 
currently are not included in the code. Thus, SUSYJ'LAVOR can be used to estimate such decay rates only 
at low tan/3 < 10. 
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However, these existing codes are restricted to the Minimal Flavor Violation scenario. Thus, 
to the best of the authors' knowledge, SUSY_FLAVOR is the first program which can simul- 
taneously calculate the set of rare decays listed in Table [1] without any (apart from the 
i?-parity conservation) restrictions on the choice of MSSM parameters. Other publicly 
available codes that are relevant to SUSY_FLAVOR (which can e.g. calculate the MSSM soft 
parameters used as input to SUSY_FLAVOR, or for the same set of input parameters cal- 
culate non-FCNC related observables) are FeynHiggs [31], SoftSUSY [32], SuSpect [33] . 
SPheno [34], MicrOMEGAs [35], DarkSUSY [36] and NMHDECAY [37]. 

In summary, the basic features of SUSY_FLAVOR are: 

• The program utilizes the most general i?-parity conserving Lagrangian for the MSSM. 
In addition to standard soft breaking terms, it can even accommodate additional non- 
holomorphic terms, such as 

A'J'Hf^QlD' + Xl'HrQl + H.c. , (1) 

that, for example, do not appear in the minimal supergravity scenario but are present 
in the most general softly broken supersymmetric effective Lagrangian [38] . 

• There is no limit on the size of flavor violating parameters because the calculation 
does not rely on the MIA expansion. Complex "mass insertions" of the form 

(/, J denote quark flavors, X, Y denote superfield chirality, and Q indicates either 
the up or down quark superfield sector, similarly for slepton superfields) are taken 
as inputs, but they only serve to conveniently parametrize the sfermion mass matri- 
ces. SUSY_FLAVOR numerically calculates the exact tree-level spectrum and mixing 
matrices, which are later used in loop calculations. 

• As an intermediate step, parton-level form factors for quark and lepton 2-, 3- and 
4-point Green functions are calculated. They are later dressed in hadronic matrix 
elements (see Table |3] in Section [3]) to obtain predictions for the physical quantities 
listed in Table [H The set of Green's functions computed by SUSY_FLAVOR as inter- 
mediate "building blocks" is quite universal and can be used by other authors to 
calculate other processes. 

• The program runs fairly quickly. On a Mac PowerBook G4 with GNU FORTRAN 
g77 it returns the output for a single parameter set run within a second. 

We note that the current SUSY_FLAVOR version does not resum higher-order corrections in 
the limit of large tan /3. Such corrections in the MFV scenario can easily dominate the 
SM result for tan /3 > 30 — 40. However, even in that case new sources of flavor violation 
often give comparable or more significant contributions than the MFV-type tan /3-enhanced 
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corrections. Thus, for large tan /3 one should perform the resummation of leading higher- 
order terms in the presence of the non-vanishing flavor violation in the sfermion mass 
matrices. Unfortunately, such resummation is not yet fully understood, i.e., although there 
are theoretical ideas [T0 |l39] - [lTj on how to resum these contributions beyond MFV, it is 
quite difficult to implement them into SUSY_FLAVOR without losing numerical stability for 
large mass insertions. Thus, even though some parts of SUSY_FLAVOR (e.g. BB mixing and 
the B — i> decays) were originally devised to perform a large tan/3 resummation for 
the MSSM parameter choice restricted to the MFV case, in the present version this option 
has been deactivated for consistency. We hope to improve this in future versions of the 
code. The current version of SUSY_FLAVOR should be used for low to moderate values of 
tau/S (tan/3 < 30), however even in the case of significant supersymmetric fiavor violation 
it should still produce reasonably accurate results even for tan (3 beyond this range. 

The rest of the paper is organized as follows. In Section [2] we define the general structure of 
the MSSM Lagrangian following Ref. [H] to facilitate comparison of the conventions used 
in SUSY_FLAVOR with others used in the literature and to connect the variables used in the 
code with physical quantities. Section |3] describes the internal structure of SUSY_FLAVOR, 
the most important steps of calculations, and the file structure of the library. In Section |4] 
we carefully present the initialization sequence for SUSY_FLAVOR, defining input parameters 
and how they are used. Routines for calculating the FCNC and CPV observables collected 
in Table [1] are described in Section [51 We conclude with a summary of the presentation. 
Appendix [A] contains brief instructions on how to install and run the SUSY_FLAVOR package. 
In Appendices |B] and [O we provide templates for initializing SUSY_FLAVOR from within the 
program and using an external file in the SLHA2 format respectively. Both of these 
templates produce the set of test results listed in Appendix [Dl 

SUSY_FLAVOR can be downloaded from the following address: 

http : //www . f uw . edu . pi/ susy_f lavor 



2 Lagrangian, conventions and the tree level masses 
2.1 Lagrangian parameters 

SUSY_FLAVOR follows the conventions for the MSSM Lagrangian and Feynman rules for the 
most general i?-parity conserving version of the MSSM Lagrangian given in [T3]. Over 100 
Lagrangian parameters are taken as input to SUSY_FLAVOR and can be initialized indepen- 
dently. 

For completeness and for easier comparison with conventions used in other sources, we 
present here the full list of the general MSSM couplings. They can be classified by sectors 
of the theory: 
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1. Gauge sector, gi, g2, gs denote the coupling constants of gauge groups U (l)y , SU (2) l, 
SU{3)c, respectively. 

2. Superpotential and Yukawa couplings. The superpotential and the soft breaking 
sfermion couplings are written after the rotations of superfields to the super-KM basis in 
which the Yukawa couplings are diagonal and the soft parameters are redefined accordingly 
to account for accommodate these flavor rotations (see e.g. [3]). 

We do not assume the existence of the heavy right-handed neutrino/sneutrino supermul- 
tiplet and neglect related terms in the Lagrangian§| Then the most general form of the 
i?-parity conserving MSSM superpotential takes the form: 

W = fie,,HlH] + e,,Y/HlL^R' + e.^Y^HlQ'p' + e^.Y^H^Qp' . (3) 

Capital indices (/, J,K...) label matter field generations and run from 1 to 3. Lower-case 
indices (z,j, ...) are SU{2)l indices (we use ei2 = —1). SU{?>)c indices are not written 
explicitly; we assume that the Q supermultiplets are QCD triplets and the D and U super- 
multiplets are anti-triplets. At tree level, quark and lepton masses are related to Yukawa 
couplings by (note that Y/, Y"/ are negative in our convention): 

It follows that in SUSY_FLAVOR the fermion masses and the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, rather than the Yukawa couplings, are used as input 
parameters. 

3. Soft gaugino mass terms for the t/(l)y, SU{2)l and SU{?>)c gauge groups 

^MiAbAb + \m2\\\\ + ^MgA^A" + h.c. (5) 

4. Soft-breaking mass terms for the scalar fields. 

-mjj^Hl^Hl - mjj^Hf*Hf - (m^ )^-^Lf*L/ - {m\y^E^*E^ 
-{KmlK^y-^Q{*Q( - imiy-^Qi*Qi - {miy-^ D^*D-^ - {miy-^U^*U^ . (6) 

5. Trilinear scalar couplings corresponding to superpotential Yukawa terms. 

ml,e,,HlH] + e,,A\' HIl]E'' + e,,A'/HlQ]D' + e,,Ai' H^Q]u' + h.c. (7) 

6. Non-standard trilinear scalar couplings involving complex conjugated Higgs fields 
(sometimes called "non-analytic terms"). 

J^i'^H^'^LlE-^ + Xl^H^^QlD' + Xl'Hl^QlU' + h.c. (8) 



^Thc modifications to the phenomenology of the MSSM from the presence of a heavy right neutrino 
supermuhiplet are discussed in [43]. Some numerical codes concerning the problem can be obtained from 
its authors. 
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Usually these couplings are not considered as they are not generated in standard SUSY 
breaking models. However, for completeness they are included in SUSY_FLAVOR and by 
default initialized to zero. Users of SUSY_FLAVOR may decide to set them to some non- 
vanishing values in order to check their impact on rare decay phenomenology. 

In general, the mass parameters fj,, m^2) '^q; ^1,2,3, and the trilinear soft couplings may 
be complex. Global rephasing of all fermion fields of the theory and of one of the Higgs 
multiplets can render two of these parameters real We choose them to be the gluino 
mass M3 and the soft Higgs mixing term m^2- The latter choice keeps the Higgs vacuum 
expectation values (VEV) and, therefore, the parameter tan /3 real at tree level. 



2.2 Physical tree level masses and mixing angles 

Mass matrices of the MSSM particles can be written in terms of the parameters of Sec- 
tion 12.11 In SUSY_FLAVOR, following [T3] , we consistently use matrix notation for all fields, 
including neutral and charged Higgs bosons. Such a notation simplifies the expressions for 
loop calculations. In this Section we explicitly write down all mass matrices to fix our sign 
conventions relative to other choices in the literature. 

1. Higgs sector. We denote the CP-even and CP-odd neutral scalars as and A^, 
respectively, with z = 1,2. In terms of more common notation, (if°,i^2) = {H^,h^) and 
[A^, A2) = (v4°, These are related to the initial Higgs doublets by (no sum over i): 

QmHl = -^ZUl (9) 

In these formulae vi, V2 are the VEVs of the two neutral components of the Higgs doublets 
and Zr, Zh are the mixing matrices in the CP-even and CP-odd Higgs sectors, respec- 
tively. 

The mixing matrix Zr and the masses of H':^ can be obtained by diagonalizing the CP-even 
Higgs mass^ matrix: 

ZT I -'"is + - % ) z - f ''^"f " \-( ° (10) 



WW 



A\{= A^) has mass M\ = mjj^ +mjj^ + 2|/ip. SUSY_FLAVOR assumes i?^ gauge with ^ 
so the neutral Goldstone boson ^^(^ mass Mqo = Mz- 



The mixing matrices Zh, Zr are parametrized as follows: 
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sin/3 —cos/3 \ Z i ^^^'^ — sina 

cos/3 sin/3 / ' ^ \ sina cosa 
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with the angles a and /3 determined by 

tan /3 = — , 

tan 2a = tan 2/3 , 



0</3<f , 



-|<a<0. 



(12) 



Charged Higgs scalars are denoted by Hf = H^,G^ and are related to the initial Higgs 
doublet again by the matrix Zh'- 



Hi* 
HI 



Ht 



The physical charged Higgs boson has mass 
while the charged Goldstone bosons have masses Mq± = M] 



(13) 



(14) 



2. Gaugino sector. The chargino masses and mixing matrices and Z_ are defined by 
the relation 




M. 



Xi 







M 



(15) 



X2 



In SUSY_FLAVOR we choose Z_, Z^ such that both masses My. are real positive and My^ > 



The neutralino tree level masses are given by 
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Ml 







ev2 
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2cw 


2cw 







M2 


evi 
2sw 


-ev2 
2sw 






—ev\ 
2cw 


evi 
2sw 





-/i 




V 


ev2 
2cw 


—ev2 
2sy^r 







/ 



V 



\ 



(16) 



where again we use the ambiguity in the definition of the Zn matrix to choose to make all 
MyO real positive and increasingly ordered. 

3. Slepton sector. The three complex sneutrino fields have tree level masses and the 
mixing matrix Zj:, defined by: 



4 ( e\vl-vl) ^ 



+ mL] Zi, 



(Ml 



\ 



\ 



Ml j 



(17) 



The mass matrix for the six charged sleptons 



/ e\v\-vl){l-2cl,) , vlY^ 



Ml 



2 \T 



4c 



w ^ 



is diagonalized by the unitary matrix Z/,, 



(Ml 



\ 



\ 



Ml J 



(19) 



4. Squark sector. Analogously, for the up and down squarks one has: 



Ml 



U 



(20) 



( Ml^ 



\ 



\ 



J 



(21) 



D 



A 



't^ 



211. /it 



+ 



(22) 



V 



\ 



Ml J 



(23) 



Note that Zu is defined with a complex conjugate compared to the definitions of Zl and 
Zj:,. Thus all positively charged sfermion mass eigenstates are multiplied by Z^^, while 
negatively charged eigenstates are multiplied by Z^*. 



2.3 Interfacing with the Les Houches Accord 

SUSY_FLAVOR has been in development since 1996, long before the Les Houches Accord 
(SLHA) for common MSSM Lagrangian conventions was agreed. Because of that, it was not 
feasible to change the internal SUSY_FLAVOR structure as it would require careful checking 
and rewriting of thousands of lines of a complicated code. Therefore we have decided to 
keep the conventions of [H] for the internal calculations in SUSY_FLAVOR. In Table [2] we 
summarize the differences of our conventions and those of the latest extended SLHA 2 [12] . 
These differences are quite minor and translation can be done by changing few signs and/or 
transposing matrices in the soft SUSY breaking sector. Thus, for the input parameters of 
SUSY_FLAVOR we leave the choice of convention as a user-defined option. 
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SLHA 2 [12] 


Ref. [H 


Td, Te 
m^, mf 


m^, ml 
(ml)^, (ml)^, {mlf 

{Mir, {Mir 



Table 2: Comparison of SLHA [12] and Ref. [H] conventions. 

Currently SUSY_FLAVOR does not use the super-PMNS basis for the lepton and slepton 
sector; only the charged lepton Yukawa matrix (and not the neutrino mass matrix) is 
diagonalized. The super-PMNS basis can become helpful once new experiments are able 
to identify the flavor of the neutrinos produced in rare decays, but at present this is not 
experimentally feasible. 

3 Structure of the code 

Calculations in SUSY_FLAVOR take the following steps: 

1. Parameter initialization. This is the most important step for SUSY_FLAVOR users 
and is described in detail Section HI Users can adjust the basic Standard Model parameters 
according to latest experimental data and initialize all (or the chosen subset of) supersym- 
metric soft masses and couplings and Higgs sector parameters listed in Section 12.11 

2. Calculation of the physical masses and the mixing angles. After setting the 
input parameters, SUSY_FLAVOR calculates the eigenvalues of the mass matrices of all MSSM 
particles and their mixing matrices at the tree level. Diagonalization is done numerically 
without any approximations. 

3. Calculation of Wilson coefficients at the SUSY scale. Physical tree-level masses 
and mixing matrices are used to evaluate exact one-loop Wilson coefficients of the effective 
operators required for a given process. Again, the formulae used in the code are exact, i.e. 
do not rely on any approximations, such as the MIA expansion. In the current version, 
SUSY_FLAVOR calculates Wilson coefficients generated by the diagrams listed in Table [31 All 
Wilson coefficients are calculated at the high energy scale, assumed to be the average mass 
of SUSY particles contributing to a given process or the top quark scale. 

It is important to stress that SUSY_FLAVOR accepts fermion generation indices and Higgs 
boson indices as input parameters. Thus in Table|3](i and u, I and v denote quarks or leptons 
of any generation and, similarly, and denote any type of the neutral Higgs bosons. 
Hence, the actual number of amplitudes which can be calculated using combinations of these 
form factors is much larger than needed for the rare decay rates currently implemented fully 
in SUSY_FLAVOR. We plan to add new processes in future releases of our library. 
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Box 


Penguin 


Self energy 


dddd 


Zdd, 'ydd, gdd 


(i- quark 


uuuu 


Hfdd, A^dd 


u- quark 


ddll 


H^uu, A^uu 




ddvp 







Table 3: One loop parton level diagrams implemented in SUSY_FLAVOR. 



4. Strong corrections. In its final step SUSY_FLAVOR performs (when necessary) the QCD 
evolution of Wilson coefficients from the high energy (SUSY or top quark mass) scale to 
the low energy scale appropriate for a given rare decay, calculates the relevant hadronic ma- 
trix elements, and outputs predictions for physical quantities. The formulae for QCD and 
hadronic corrections are primarily based on calculations performed in the SM and supple- 
mented, when necessary, with contributions from non-standard operators which usually are 
neglected in the SM, because they are suppressed by powers of the light quark Yukawa cou- 
plings. This part of SUSY_FLAVOR is based on analyses published by other authors, whereas 
points 1-3 are implemented using our own calculations. The accuracy of strong corrections 
differ from process to process, from negligible or small (leptonic EDM, "gold-plated" de- 
cay modes K — ttuu [45]) to order of magnitude uncertainties (unknown long distance 
contributions to ArriK or Amo)- Even in the case of large QCD uncertainties, the result 
of the calculation performed by SUSY_FLAVOR can be of some use. Flavor violation in the 
sfermion sector can lead to huge modifications of many observables, sometimes by several 
orders of magnitude, so that comparison with experimental data can help to constrain the 
soft fiavor-violating terms even if strong corrections are not very well known. 

Below we list the files included in the SUSY_FLAVOR library with a brief description of their 
content and purpose. 



eischl 


f : 


auxiliary numerical routine - hermitian matrix diagonalization 


vegas 


f : 


auxiliary numerical routine - Vegas Monte Carlo integration 


rombint 


f : 


auxiliary numerical routine - Romberg numerical integration 


sf lav_io 


f : 


input routine for reading of the SLHA2 format; test output routines 


b_f un 


f : 


general 2-point loop functions 


db_f un 


f : 


derivatives of general 2-point loop functions 


c_f un 


f : 


general 3-point loop functions 


cd_f un 


f : 


3-, 4- and some 5-point loop functions at vanishing external momenta 


vh_def 


f : 


definitions of Higgs boson tree-level vertices 


vg_def 


f : 


definitions of gauge boson tree-level vertices 
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vf_def .f definitions of fermion tree-level vertices 

mh_init . f : initialization of MSSM parameters 

mh_diag . f : diagonalization of tree level mass matrices; outputs physical masses and mixing angles 

qcd_fun.f : auxiliary QCD calculations - running a^, running quark masses etc. 

d_self . f : (i-quark self-energy 

u_self . f : M-quark self-energy 

sff_fun0.f : form factors of the general scalar-fermion-fermion 1-loop triangle diagram 

sdd_vertO . f : CP-even Higgs-d quark-d quark 1-loop triangle diagram 

pdd_vertO . f : CP-odd Higgs-c? quark-d quark 1-loop triangle diagram 

suu_vertO . f : CP-even Higgs-M quark-w quark 1-loop triangle diagram 

puu_vertO . f : CP-odd Higgs-M quark-w quark 1-loop triangle diagram 

zdd_vertO . f : Z boson-c? quark-rf quark 1-loop triangle diagram 

ddg_fun.f : form factors for the general gauge boson-fermion-fermion 1-loop triangle diagram 

dd_gluon . f : d quark-rf quark-gluon 1-loop triangle diagram 

dd_gainma . f : d quark-c? quark-photon 1-loop triangle diagram 

bsg_nl .f : formulae for Br (5 — )■ Xs7), including QCD corrections 

dd_ll . f : d quark-c? quark- lepton-lepton 1-loop box diagram 

dd_vv . f : d quark-c? quark-neutrino-neutrino 1-loop box diagram 

phen_2q.f : formulae for Br{Kl -> tt^uu), Br(i^+ vr+Pi/) and Bi{Bs(^d) -> including QCD 

corrections and hadronic matrix elements 

dd_mix . f : 4-d quark 1-loop box diagram 

uu_mix . f : 4-u quark 1-loop box diagram 

phen_4q.f : formulae for the meson mixing observables: Atuk, ^k, Am/), Am^^^^j including QCD 
corrections and hadronic matrix elements 

edin_l . f : lepton electric dipole moment 

cdin_d . f : (i-quark chromoelectric dipole moment 

cdm_u . f : n-quark chromoelectric dipole moment 

cdni_g . f : gluon chromoelectric dipole moment 

edin_d . f : d-quark electric dipole moment 
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edm_u . f : ti-quark electric dipole moment 

edm_n . f : n-quark electric dipole moment 

All the 2-, 3- and 4-point Green functions are calculated for vanishing external momenta. 
As mentioned before, by "u quark" and "c? quark" we mean all generations of quarks. 

In addition to the files listed above, the library contains the master driver file susy_f lavor . f 
which illustrates the proper initialization sequence for SUSY_FLAVOR parameters and pro- 
duces a set of test results for the implemented observables. 

4 Parameter initialization in SUSY_FLAVOR 

We now list the input parameters used by SUSY_FLAVOR. These are not always directly the 
MSSM Lagrangian parameters given in Section 12.11 - for example, instead of using the /i 
parameter and the soft Higgs masses m^.,m^2) is customary to use tan/3 and the CP- 
odd Higgs mass Ma to parametrize MSSM Higgs sector. In its first step, SUSY_FLAVOR 
restores the Lagrangian parameters of Section l2.1l for the given set of more human-friendly 
input parameters. Then, the remaining routines use the "raw" Lagrangian parameters — if 
necessary they can also be directly modified by (experienced!) users. 

In the rest of this section we describe step-by-step the basic initialization routines used by 
SUSY_FLAVOR, their arguments and, when necessary, the FORTRAN common blocks storing 
the most important data (other common blocks serve for the internal purposes and usually 
do not need to be accessed by users). 

By default, SUSY_FLAVOR uses the following implicit type declaration in all routines: 
implicit double precision (a-h,o-z) 

so that all variables in SUSY_FLAVOR with the names starting from a to h and from o to z are 
automatically defined as double precision and those with names starting from i to n are 
of integer type. In what follows we explicitly indicate variables that do not obey this rule. 
Such variables are always listed in explicit type statements inside the procedures. Complex 
parameters mentioned in this article are declared in SUSY_FLAVOR as double complex type. 
Mass parameters are always given in GeV. 

SUSY_FLAVOR provides two ways of initializing input parameters. As the first option, they 
can be read from the file susy_f lavor . in. The structure of this file follows the SLHA2 con- 
vention |12], with some extensions which we describe in Section W7[\ Initializing parameters 
in the input file is simple, it is done by a call to single subroutine sf lav_input and does not 
require detailed knowledge of the program internal structure. This option is particularly 
convenient for testing a single parameter set but can be cumbersome for scans over the 
MSSM parameter space. Therefore, as a second option, SUSY_FLAVOR also provides a set of 
routines designed to initialize parameters defined in the program, which can easily be used 
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to prepare programs that scan over large parameter sets. As described in Section I4.2[ these 
routines require more care in use, as they should be initialized in the proper order, i.e. the 
gauge sector first, then the fermion sector, Higgs sector, and SUSY sectors at the end (the 
initialization sequences for the gaugino, slepton and squark sectors are independent). 

An example of a full initialization sequence for SUSY_FLAVOR, illustrating both options 
mentioned above, is presented Appendix [Bl The sample input file susy_f lavor . in is given 
in Appendix O Test output generated for parameters used in Appendices [B] and [C] is 
enclosed in Appendix 



4.1 Parameter initialization from the input file 

Input parameters for SUSY_FLAVOR can be set by editing appropriate entries of the file 
susy_f lavor . in and subsequently calling the subroutine sflav_input, which reads the 
input file, stores the the MSSM Lagrangian parameters in FORTRAN common blocks and 
calculates tree-level physical masses and mixing matrices. After calling sf lav_input, all 
physical observable described in Section |5] can be calculated. 

The input file susy_f lavor . in is written in the SLHA2 format, with some extensions which 
we list below (for an example of a complete input file see Appendix [C|) . 

1. We define a non-standard Block SOFTINP. Currently it contains two control variables, 
iconv and input_type. These serve to choose input conventions in the sfermion sector (in 
other sectors SLHA2 and Ref. |T3] agree). 

Variable value Sfermion sector parametrization 

iconv = 1 MSSM parameters defined in SLHA 2 conventions. 

iconv = 2 MSSM parameters defined in conventions of Ref. 

input_type = 1 sfermion diagonal trilinear mixing terms given as dimensionless pa- 
rameters; all off-diagonal soft terms are given as dimensionless mass 
insertions — see comments below on the data blocks defining the 
sfermion soft terms. 

input_type = 2 sfermion soft terms given as absolute values of dimension mass^. 

2. SUSY_FLAVOR uses the W boson mass as a basic parameter rather than Fermi constant 
Gp. Therefore, in susy_f lavor . in we use entry 30 of Block SMINPUTS (not used in the 
standard SLHA2) to define M^. 

3. We allow complex values for /i and two of the gaugino masses — chosen to be the U{1) and 
SU (2) mass terms Mi and M2. Their real and imaginary parts are defined in blocks EXTPAR 
and IMEXTPAR. We use tan /3 and the CP-odd Higgs mass M4 as the input parameters for 
the Higgs sector. 

4. Following the SLHA2 convention, we only define the upper triangle of each of the 
hermitian sfermion soft mass matrices in the MSL2IN, MSE2IN, MSq21N, MSD21N, MSU2IN 
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and IMMSL2IN, IMMSE2IN, IMMSq2IN, IMMSD2IN, IMMSU2IN blocks. It is obligatory to 
define all entries, both diagonal and upper off-diagonal, since SUSY_FLAVOR does not read 
diagonal sfermion masses from the EXTPAR block. The iconv parameter defined in the 
SOFTINP block determines if sfermion parameters are given in SLHA2 or Ref. [2] conven- 
tions (see Table [2]). Finally, the input _type parameter in the SOFTINP block defines the 
format of the off-diagonal mass terms. If input_type = 1, the off-diagonal entries given in 
susy_f lavor . in are assumed to be dimensionless mass insertions and the flavor violating 
sfermion mass terms are calculated as 

iOu = (On = ^'x'\Jimjc)n{ml)jj , (24) 

where X = L,E, Q, U, D and /, J enumerate superpartners of the mass-eigenstate quarks. 

5. The blocks TEIN , TDIN , TUIN and IMTEIN , IMTDIN , IMTUIN define the trilinear sfermion 
mixing matrices which are generally non-hermitian. One is required to define all entries. 
As for the soft mass terms, the iconv parameter chooses the input convention, SLIIA2 
or Ref. [H]. For the trilinear mixing, the parameter input_type defines the format and 
dimension of both the diagonal and off-diagonal terms. If input_type = 1, then all rel- 
evant susy_f lavor . in entries are treated as dimensionless numbers and expanded to full 
trilinear mixing matrices using eqs. ( 125|26|) . For the diagonal LR terms, SUSY_FLAVOR uses 
the formulae 




where af , a^, are the diagonal trilinear mixing terms read from the input file. 
For the off-diagonal LR terms, SUSY_FLAVOR uses 

= ^hiR^Hi^DiiiOjJ ■ (26) 

Note that in eqs. f l2^|2^ for simplicity we use {TnQ)jj as the diagonal mass scale for both up 
and down left squark fields (in general related by the CKM rotation, see eqs. ( 120)1221) ). 

4.2 Parameter initialization inside the program 

SUSY_FLAVOR input parameters can also be initialized directly inside the driver program 
using the set of routines described below. Before the proper initialization sequence, the 
user can set the iconv variable value to choose the input convention: 
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common/ sf _cont/eps , indx (3,3), iconv 

iconv=l SLHA2 [12] input conventions 

iconv=2 [13] input conventions 

After choosing the input conventions, one should subsequently initialize the gauge, matter 
fermion, Higgs, SUSY fermion and sfermion sectors, using the procedures described in detail 
in the following sections. 



4.2.1 Gauge sector 

As input, SUSY_FLAVOR takes the gauge boson masses {Mw,Mz) and the gauge coupling 
constants (electromagnetic and strong) at the Mz scale. They are initialized by: 

Routine and arguments Purpose and MSSM parameters 

vpar_update(zm,wm,alpha_em) Sets electromagnetic sector parameters 

zm Mz, Z boson mass 

wm Mw, W boson mass 

alpha_em aem{Mz), QED coupling at Mz scale 

laiii_f it (alpha_s) Sets as{Mz) and Aqqd for 4-6 flavors at the NNLO 

level 

laiii_f it_nlo (alpha_s) Sets as{Mz) and Aqcd for 4-6 flavors at the NLO 

level 

alpha_s as{Mz), strong coupling at Mz scale 



4.2.2 Matter fermion sector 

SUSY_FLAVOR assumes that neutrinos are massless. Pole masses of the charged leptons are 
initialized in the file mh_init.f in block data init_phys. They are stored in the em 
array in common/fmass/em(3) ,um(3) ,dm(3) and can be directly modified there. Their de- 
fault values are: 
Lepton mass Value 

me em(l) = 0.000511 

em(2) = 0.105659 
ttIt- em(3) = 1.777 

In the quark sector the most important input parameters are the running top and bottom 
masses at a given renormalization scale and the CKM matrix angles and phase. They can 
be set by: 
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Routine and arguments Purpose and MSSM parameters 

init_f erniion_sector (tm,tscale ,bm,bscale) Sets running top and bottom quark mass 
tm,tscale mt{fit), running MS top quark mass 

biii,bscale mi,{fii,), running MS bottom quark mass 

ckm_init(sl2,s23,sl3,delta) Initialization of the CKM matrix 

sl2,s23,sl3 sin ^12, sin ^23; sin ^13, sines of the CKM 

angles 

delta S, the CKM phase in radians 

The light quark masses are also initialized in the block data init_phys of the file mh_init . f 
and stored in common/fmass_high/umu(3) ,uml(3) ,amuu(3) ,dmu(3) ,dml(3) ,amud(3). The 
arrays uml(dml) contain up(down) quark masses at the scale amuu(amud), respectively. 
Their default values are: 



Running quark mass 


Mass value 


Mass scale 






dml(l) = 0.007 


amud(l) = 


2 




dml(2) = 0.11 


amud(2) = 


2 




dml(3) = 4.17 


amud(3) = 


4.17 




uml(l) = 0.004 


amuu(l) = 


2 




uml(2) = 1.279 


amuu(2) = 


1.279 




unil(3) = 163.5 


amuu(3) = 


163.5 



The variables of the arrays uml , amuu , dml , amud can be directly accessed and modified 
if necessary. However, for consistency, after such modifications the user should call the 
routine init_run_qmass which calculates running quark masses at the high rrit scale (stored 
in common/fmass_high/ in the arrays umu,dmu and in common/fmass/ in the arrays um,dm) 
for later use in the running Yukawa couplings and in SUSY loop calculations. 



4.2.3 Higgs sector 



Following the common convention, we take the Higgs mixing parameter the CP-odd 
Higgs boson mass M^, and the ratio of vacuum expectation values tan /3 = V2/V1 as the 
input parameters. Other Higgs sector parameters listed in Section 12.11 can be expressed 
as: 

ml, = ^(Ml-2|/ip-(Ml + M|)cos2/3), 
= i(Ml - 2|/i|2 + (Ml + M|) cos 2/3), 

= (27) 

^2 2sin2/3 ^ ^ 

The MSSM Higgs sector at the tree level can be effectively parametrized in terms of just Ma 
and tan f3, but the fi parameter is necessary for the chargino and neutralino sectors. Here it 
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is used to calculate the original Higgs soft mass parameters m|^^ and fnjj^ for completeness 
and future applications; they currently have no further use. 



Routine and arguments 

init_higgs_sector (pm , tb , amu , ierr) 
pm 
tb 
amu 
ierr 

init_yukawa 



Purpose and MSSM parameters 

Higgs sector initialization 

CP-odd Higgs mass 

Ratio of Higgs VEVs, tan /3 = ^ 

Higgs mixing parameter fi (complex) 

output error code: ierr 7^ if Higgs sector 

initialization failed 

Initialization of the running Yukawa cou- 
plings Yi,Yu,Yd for all generations (at the 
same scale as the running quark masses) 



4.2.4 Supersymmetric fermion sector 



Initialization is done by the routine: 
Routine and arguments 

init_ino_sector (gml , gm2 , gm3 , amu , tb , ierr) 
gml , gm2 
gm3 
tb 
amu 
ierr 



Purpose and MSSM parameters 

gaugino sector initialization 

U{1), SU{2) gaugino masses (complex) 

SU{3) gaugino mass 

tan /3 = ^, the ratio of Higgs VEVs 

the Higgs mixing parameter fi (complex) 

output warning code: ierr 7^ if a char- 

gino or a neutralino is lighter than Mz/2 



If one sets Mi = in the call to init_ino_sector then the GUT-derived relation Mi 
I tan^ 6*1^^2 is used in the gaugino mass calculations. 



4.2.5 Sfermion sector 

This is the most complicated MSSM sector; it contains a large number of free param- 
eters. SUSY_FLAVOR supplies two subroutines for the sfermion parameters initialization, 
init_slepton_sector and init_squark_sector. They accept as input only dimension- 
less mass insertions and dimensionless diagonal trilinear soft mixing terms, expanded in 
SUSY_FLAVOR to entries of the soft mass matrices as defined by eqs. fl24|25|26l) (this is only 
a particular choice of parametrization and does not lead to any loss of generality). The 
sfermion initialization routines have the following arguments: 
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subroutine init_slepton_sector (sll , sir , asl , ierr , slmi_l , slmi_r) 



Argument 

sll 
sir 
asl 



lerr 
slmi_l 



slmi_r 



slmi_lr 



MSSM parameters 

Array of the diagonal left slepton masses {m\)ii = sll(l)^, J = 1...3 
Array of the diagonal right slepton masses (m|;)// = slr(l)^, / = 1 . . . 3 
Array of the dimensionless diagonal slepton trilinear mixing terms aj = 
asl(l),/ = 1...3 (complex parameters). 

output error code: ierr 7^ if slepton sector initialization failed 

Array of the off-diagonal left slepton mass insertions 6}^ = slnii_l(l), 

S]^ = slnii_l(2), 6}^ = slmi_l(3) (complex parameters); remaining LL 

mass insertions are initialized via hermitian conjugation 

Array of the off-diagonal right slepton mass insertions 6^ = slmi_r(l), 

5|f = slmi_r(2), = slini_r(3) (complex parameters); remaining RR 

mass insertions are initialized via hermitian conjugation 

Matrix with off-diagonal slepton trilinear LR mass insertions (Jf*^^ = 

slnii_lr(l, J), J, J = 1 ... 3 (complex parameters) 



subroutine init_squark_sector (sql , squ , sqd , asu , asd , ierr , 
sqini_l , sunii_r , sdnii_r , sumi_lr , sdnii_lr) 

MSSM parameters 

Array of the diagonal left squark masses {m?Q)ii = sql(l) 



Argument 

sql 
squ 
sqd 

asu 

asd 

ierr 
sqmi_l 

sumi_r 

sdmi_r 

sumi_lr 
sdmi_lr 



Array of the diagonal right up-squark masses (m^)// = squ(l)^, 



= 1...3 
1 = 1.. 
sqd(l)', / 



Array of the diagonal right down-squark masses 
1...3 

Array of the dimensionless diagonal soft LR up-squark mixing terms = 
asu(l), / = 1 . . . 3 (complex parameters) 

Array of the dimensionless diagonal soft LR down-squark mixing terms 

= asd(l), / = 1 . . . 3 (complex parameters) 
output error code: ierr 7^ if squark sector initialization failed 
Array of the off-diagonal left squark mass insertions = sqini_l(l), 



= sqnii_l(2 



6^ = sqini_l(3) (complex parameters); remaining QLL 
mass insertions are initialized via hermitian conjugation 
Array of the off-diagonal right up-squark mass insertions Sj^ = sunii_r(l), 
6ff = sumi_r(2), 6jf = sunii_r(3) (complex parameters); remaining URR 
mass insertions are initialized via hermitian conjugation 
Array of the off-diagonal right down-squark mass insertions 6}^ = 
sdnii_r(l), 6]^ = sdini_r(2), 6]^ = sdnii_r(3) (complex parameters); re- 
maining DRR mass insertions are initialized via hermitian conjugation 
Matrix with off-diagonal up-squark trilinear LR mass insertions ^[/lji = 
sumi_lr(l, j), /, J = 1 . . . 3 (complex parameters) 

Matrix with off-diagonal down-squark trilinear LR mass insertions Sj^j^p^ = 
sdmi_lr(l, J), J, J = 1 ... 3 (complex parameters) 
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If necessary, experienced SUSY_FLAVOR users can directly modify the soft breaking sfermion 
parameters stored in common blocks /msoft/ and /soft/ (see Table|4]in Section I4l3l) . One 
must subsequently call the routines sldiag, sqdiag (see file mh_diag.f) to recalculate 
the tree-level sfermion masses and mixing matrices. This may, however, require a deeper 
understanding of the SUSY_FLAVOR initialization sequence and its data structure. 

4.3 Tree-level physical masses and mixing angles 

After performing the full initialization sequence in SUSY_FLAVOR, all the MSSM Lagrangian 
parameters listed in Section 12.11 physical tree-level particle masses (with the exception 
of the running quark masses), and mixing matrices are calculated and stored in common 
blocks. If necessary, they can be directly accessed and modified. Note, however, that after 
any modifications of the Lagrangian parameters, relevant procedures calculating physical 
masses and mixing angles have to called again. In Table H] we list the important blocks 
storing MSSM parameters. Common blocks containing masses and mixing angles are listed 
in Table [51 

All parameters, tree-level masses, and mixing angles can be printed for test purposes, e.g. 
by calling the subroutines print_MSSM_par and print_MSSM_masses. 

5 List of processes 

In this section we list the set observables whose computation is fully implemented in 
SUSY_FLAVOR vl.O. For all of them, SUSY_FLAVOR takes into account one-loop supersym- 
metric contributions. QCD corrections and hadronic matrix elements are extracted from 
the papers of various authors, mostly from analyses done in the Standard Model. They are 
assumed to work reasonably well also in the MSSM since supersymmetric strong corrections 
from gluino and squarks are suppressed by large masses of these particles. 

In most cases, QCD and hadronic corrections are known at the level of few to tens %, while 
variations of supersymmetric flavor and CP violating parameters can change observables by 
orders-of-magnitude. Thus, as long as the MSSM parameters are not measured very pre- 
cisely, the current implementation of strong corrections is sufficient for analyses performed 
in the framework of the general MSSM. 

Calculations of the hadronic matrix elements are particularly difficult as they have to be 
performed, at least partially, in the regime of strongly coupled QCD. Results of such cal- 
culations can differ significantly depending on the methods used and thus carry significant 
theoretical uncertainties. Therefore, in SUSY_FLAVOR, quantities which requires hadronic 
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Common block and variables 



Lagrangian parameters 



common/vpar/ st , ct , st2 , ct2 , set , sct2, e , e2 , alpha, wm,wm2 ,zm,zm2 , pi , sq2 



St , ct , st2 , ct2 , set , sct2 

e , e2 , alpha 
wm , wm2 , zm , zm2 
pi , sq2 

common/hpar/hml , hm2 , hml2 , hmu 
hml ,hm2 
hml2 
hmu 

common/vev/vl ,v2 
vl,v2 

common/yukawa/yl(3) ,yu(3) ,yd(3) 
yl(3) 
yu(3) 

yd(3) 

common/ gmass/gm3 , gm2 , gml 
gml,gm2 
gm3 



Weinberg angle functions, respectively sw, cw, 



electric charge powers at Mz scale: e, e^, ctem 
gauge boson masses: Mw, M^i, Mz, M| 
numerical constants, n and ^/2 



soft Higgs masses m|^^ , mfj^ 

soft Higgs mixing parameter 

Higgs mixing parameter /i (complex) 

Higgs vacuum expectation values vi,V2 

charged lepton Yukawa couplings Ye, Y^, Y^ 
Running MS up-quark Yukawa couplings at rrit 
scale: Fc, Yt 

Running MS down-quark Yukawa coupUngs at 
mt scale: Yu,Yc, Yt 

U{1),SU{2) gaugino masses Mi,M2 (complex) 
SU (3) gaugino mass M3 



common/msoft/lms(3,3) ,rms(3,3) ,ums(3,3) ,dms(3,3) ,qms(3,3) 

lms(3,3) ,rms(3,3) hermitian slepton soft mass matrices m\, m% 

(complex) 

ums(3,3) ,dms(3,3) ,qms(3,3) hermitian squark soft mass matrices mfj, mj^, 

ttlq (complex) 

common/soft/ls(3,3) ,ks(3,3) ,ds(3,3) ,es(3,3) ,us(3,3) ,ws(3,3) 

ls(3,3) ,ds(3,3) ,us(3,3) trilinear soft LR mixing matrices Ai, Ad, Au 

(complex) 

ks(3,3),es(3,3),ws(3,3) trilinear "non-holomorphic" soft mixing matri- 

ces A'l, A'^, A'^ (complex) 

Table 4: Common blocks storing the MSSM Lagrangian parameters. We omit flavor indices 
in the fermion and sfermion sectors. 
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Common block and variables Masses and mixing matrices 

common/f mass/em (3) ,um(3) ,dni(3) 

em(3) Charged lepton pole masses rrie, m^, rrir 

um(3) Running MS up-quark masses at the rrit scale: 

dm (3) Running MS down-quark masses at the mt 

scale: mujmcmt 

common/hmass/cm(2) ,rm(2) ,pm(2) ,zr(2,2) ,zh(2,2) 

rni(2) neutral CP-even Higgs masses rm(l) = Mh, 

rm(2) = Mh 

pni(2) neutral CP-odd Higgs mass pm(l) and Gold- 

stone mass pm(2) 

cm (2) charged Higgs mass cm(l) and charged Gold- 

stone mass cm (2) 
zr(2,2) CP-even Higgs mixing matrix Zji 

zh(2,2) CP-odd and charged Higgs mixing matrix Zh 

common/charg/f cm(2) ,zpos(2,2) ,zneg(2,2) 

fern (2) chargino masses M +, i = 1,2 

zpos(2,2) ,zneg(2,2) chargino mixing matrices Z^,Z_ (complex) 

common/neut/f nm(4) ,zn(4,4) 

f nm(4) neutralino masses M o, i = 1 . . .4 

zn(4,4) neutralino mixing matrix Z^v (complex) 

common/slmass/vm(3) ,slm(6) ,zv(3,3) ,zl(6,6) 

vm(3) sneutrino masses M^^, / = 1 . . . 3 

slm(6) charged slepton masses ML.,i = 1 ... 6 

zv(3,3) sneutrino mixing matrix Z^, (complex) 

zl(6,6) charged slepton mixing matrix Z^ (complex) 

common/sqmass/sum(6) ,sdm(6) ,zu(6,6) ,zd(6,6) 

sum (6) up-squark masses Mu-,i = 1 ... 6 

sdm(6) down-squark masses M^^, -i = 1 ... 6 

zu(6,6) up-squark mixing matrix Zu (complex) 

zd(6,6) down-squark mixing matrix Zd (complex) 

Table 5: Common blocks storing physical particle masses and mixing matrices. 
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matrix element estimates and other QCD related quantities are treated as external param- 
eters. They are initialized to the default values listed below for each observable and can be 
directly modified by users by changing the relevant variables in the common blocks where 
they are stored. Currently most of the hadronic (and related) input parameters used in 
SUSY_FLAVOR are taken from the Table 3 of Ref. [16]. 



5.1 Electric Dipole Moments of charged leptons 

Lepton EDMs are defined as the coefficient dii in the effective Hamiltonian for the flavor- 
diagonal lepton-lepton-photon interaction: 

= '-^Pa.^l^l'F'^'' , (28) 

where / = 1,2,3 is the generation index of the lepton as usual. In SUSY_FLAVOR lepton 
EDM is calculated by: 

Routine: double precision function edm_l(l) 

Input: / = 1,2,3 for e, /i, r respectively 

Output: EDM for the charged lepton specifled by I 

QCD related factors: none, QCD corrections are small and not included 

Details of calculations: Ref. |1] 



5.2 Neutron Electric Dipole Moment 

The neutron EDM can be approximated by the sum of the electric dipole moments of 
the constituent d and u quarks plus contributions of the chromoelectric dipole moments 
(CDM) of quarks and gluons. The EDMs of the individual quarks are deflned analogously 
to eq. fl25]) . The CDM Cg of quark q is deflned as: 

He = -^qa,,^,T-qG^^''\ (29) 

The gluonic dipole moment Cg is deflned as: 

T^g = ~~^fabcG°^pG^J'G\^e^''^". (30) 

The exact calculation of the neutron EDM requires knowledge of its hadronic wave function. 
SUSY_FLAVOR uses the "naive" chiral quark model approximation |47] : 

En = |(4rf. - rf„) + ^(4q - c„) + ^^c, (31) 
3 Aw An 

where rji and are the QCD correction factors [H] and the chiral symmetry breaking 
scale [17], respectively. Various models give signiflcantly different rji factors. As a result, 
even the sign of the neutron EDM is not certain. Thus the SUSY_FLAVOR result should be 
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treated as an order of magnitude estimate only. The calculations are performed by calling 

Routine double precision function edm_n() 

Input none 

Output neutron EDM 
QCD related factors: 

common/edm_qcd/eta_e , eta_c , eta_g , alamx 

r]e eta_e = 1.53 

rjg eta_c = 3.4 

% eta_g = 3.4 

Ax alamx = 1.18 

Details of calculations: Ref. |1] 

5.3 — > TT^uu and — > TT^uiy decay rates 

The relevant part of the effective Hamiltonian generated by the top quark and SUSY particle 
exchanges can be written as 

G (y. 

'Hes = ^- — [^Lisd)v-Aii^ii^i)v-A + XR{sd)v+A{i^ii^i)v-A\- (32) 
The branching ratios for the K — j- ttuu decays are then given by 



(33) 



Br(A-o^.°..) = .,( ^"^(y^^) )^ (34) 

where k [H], A (one of the Wolfenstein parameters [SD]), and the NLO charm quark con- 
tribution Pc can be modified by SUSY_FLAVOR users (note that k and Pc depend 
on Vusi and a^) Branching ratio calculations are performed by calling 

Routine subroutine k_pivv(br_kO,br_kp) 

Input none 
Output br_kO = Br{Kl tx^vv) 

br_kp = Br{K^ tt+z/i/) 

QCD related factors 

common/kpi vv/ akO , del_akO , akp , del_akp , pc , del_pc , alam 

kl ± Akl akO = 2.231 ■ 10^^°, del_akO = 0.013 ■ 10^^° 

K+ ± Afi;+ akp = 5.173 ■ 10-i\ del_akp = 0.025 ■ lO^^^ 

Pc ± APc pc = 0.41, del_pc = 0.03 

A alam = 0.225 

Details of calculations: Ref. [H] 
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5.4 53 l^+l^- and 5j l^+l^- decay rates 

The general expression for these branching ratios are rather comphcated and can be found 
in [12]. For most users it is sufficient to know that, in addition to the MSSM parameters, 
the dilepton B decays depend on the B meson masses and the hadronic matrix elements of 
the down quark vector and scalar currents: 

(0|67^PL(R)s|S,(rf)(p)) = ) (35) 

i Ml fe 

where is the momentum of the decaying i?s(d)-meson of mass Mb^^^^s^ ■ The 5° — l^^l"^' 
and i?^ — 7- l^^V-'~ decay branching ratios are calculated by: 

Routine double precision function b_ll (K,L, I , J) 

Input /, J = 1, 2, 3 - outgoing leptons generation indices 

K,L - generation indices of the valence quarks of the 
B^ meson: setting {K,L) = (3, 1), (1, 3), (3, 2) and (2,3) 
chooses respectively B^, B^, B^ and -B° decay 

Output Branching ratios of the decay defined by K, L, I, J 

QCD related factors 

common/meson_data/dmk,amk,epsk,fk,dmd,aind,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

Mb^ ajnb(l) = 5.2794 

Mb, amb(2) = 5.368 

fB, fb(l) = 0.2 

/b. fb(2) = 0.245 

Details of calculations: Ref. [12] 



5.5 K^K' meson mixing parameters 

SUSY_FLAVOR calculates two parameters measuring the amount of CP-violation in neutral 
K meson oscillations: Ek and the — mass difference AMk- 

AMk = 2Ke(i?°|i/,^^=V^°) , (37) 
exp(z7r/4) 



£k 



V2AAI 



-Qm{K''\H^/=^\K') . (3J 



K 



QCD dependent corrections are known with reasonable accuracy for the Ek parameter. The 
long distance contributions to AM^ are large and difficult to control. Thus the result given 
by SUSY_FLAVOR for AMk should be treated as an order of magnitude estimate only. 
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Apart from the MSSM parameters, the calculation of the K^K^ meson mixing requires 
knowledge of the meson masses and of the hadronic matrix elements of the following set of 
four-quark operators: 




{qil.PLqiWprPLqi), 
{qiPLqiW0PRqi), 



{qiPLqiWpPLqi), 
(qia^^PLq^Wpa^'PLqi) 



(39) 



where a, /3 are color indices, for the K^K^ mixing one should choose flavor indices 1 = 2 
and J = 1. The matrix elements can be written as: 



where Fk is the i^'- meson decay constant. By default, SUSY_FLAVOR uses the Bf values at 
the scale fi = 2 GeV given in [53] using the NDR renormalization scheme (quark masses at 
the scale 2 GeV are stored in common/fmass_high/, see Section r4.2.2p . 

In addition to the hadronic matrix elements, QCD corrections depend also on the "r;" factors 
describing the evolution of the relevant Wilson coefficients from the high to low energy 
scale. These factors are automatically calculated at NLO by SUSY_FLAVOR. For the SM 
contribution to the Wilson coefficient of the Q^^^ operator a separate careful calculation 
of the evolution factors has been performed [5ll|55]. Therefore SUSY_FLAVOR treats this 
contribution separately, setting -B^jf" and the rjsM factor to default values given in 
(see [53] for a very detailed discussion of the structure of the QCD corrections in B^B^ and 
K^K^ systems, including their renormalization scheme dependence and calculations of the 
evolution factors "r/" implemented in SUSY_FLAVOR). 

The kaon mass difference AM^ and the Sk parameter measuring the amount of CP viola- 
tion in K^K^ mixing are calculated by 




(40) 
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Routine subroutine dd_kaon(eps_k,delta_mk) 

Input none 

Output eps_k = ek parameter 

deltajnk = AM^ mass difference 

QCD related factors: 

common/meson_data/dmk,amk,epsk,fk,dmd,aind,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 



Mk 


amk = 0.497672 


Measured AM^^ 


dmk = 3.49 ■ 10" 


Measured 


epsk = 2.26 ■ 10 


fx 


fk = 0.1598 


common/bx_4q/bk(5) ,bd(5),bb(2,5) , ainu_k , amu_d , amu_b 




bk(l) = 0.61 




bk(2) = 0.76 


Bl'^Hf^K) 


bk(3) = 0.51 


BY'M 


bk(4) = 0.96 


B^^'M 


bk(5) = 1.30 


Renormalization scale hk 


amu_k = 2 


common/sm_4q/eta_cc , eta_ct , eta_tt , eta_b , bk_sni , bd_sm 


rVLL 
^SM 


bk_sni = 0.724 


Vcc 


eta_cc = 1.44 


Vet 


eta_ct = 0.47 


Vtt 


eta_tt = 0.57 


Details of calculations: 


Ref. mm 



5.6 D D meson mass difference 



Calculations of the mass difference Am^) of the neutral D mesons have large theoretical 
uncertainties due to unknown long-distance strong corrections. Thus, as in the case of 
ArriK, the SUSY_FLAVOR result for Amo should be treated as an order of magnitude estimate 
only. 

The structure of strong corrections is analogous to those in the K meson system. However, 
in this case hadronic matrix elements and QCD evolution calculations available in the lit- 
erature are much less refined. SUSY_FLAVOR uses the NLO evolution for the "77" factors and 
sets, by default, all the relevant hadronic matrix elements Bi = 1, i.e. it uses the "vacuum 
saturation" approximation (this can be changed easily when new results become available). 
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Routine subroutine uu_bmeson(delta_md) 

Input none 

Output deltajnd = AMd mass difference 

QCD related factors: 

common/meson_data/dink,amk,epsk,fk,dmd,aind,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

Md amd= 1.8645 

Measured AM^^ dmd = 4.61 ■ IQ-^^ 

fo f d = 0.165 

common/bx_4q/bk(5) ,bd(5),bb(2,5) , ainu_k , amu_d , amu_b 

BY'^'^ifiD) bd(l) = 1 

bd(2) = 1 

L(/ic) bd(3) = 1 

Bl^'ifxn) bd(4) = 1 

bd(5) = 1 

Renormalization scale fi^ amu_d = 2 

common/ sm_4q/ eta_cc , eta_ct , eta_tt , eta_b , bk_sm , bd_sm , bb_sm (2) 

Details of calculations: Performed by authors, unpublished 



5.7 B'^B^ and B^B'^ mass differences 

Mixing and CP violation phenomena are also observed in the neutral B meson systems. 
In particular, the mass differences in the B^B^ and B^B^ oscillations have been mea- 
sured, 

AM^,^^, =2|(5,%|i/,^/=2|5,%)|. (41) 

In addition to the MSSM parameters, theoretical calculations of Atub^ and Am^^ depend, 
as for K and D oscillations, on the relevant hadronic matrix elements and QCD evolution 
factors. The formulae for B^B'^ mixing can be obtained by making the obvious replacements 
in the formulae presented in Section 15.51 Currently SUSY_FLAVOR uses the same set of Bi 
factors for both the i?^ and Bg sectors, but it leaves the possibility to distinguish between 
them in future, if necessary. For this one needs to independently initialize the arrays 
bb(l,i) {Bet nieson hadronic matrix elements) and bb(2,i) [Bg meson hadronic matrix 
elements) stored in common/bx_4q/. 

The values of the B meson masses and coupling constants are the same as those listed in 
Section WM. SUSY_FLAVOR calculates the mass differences AM^^^^^ as defined by eq. ( HTl) : 
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Routine subroutine dd_bmeson(i ,delta_mb) 

Input i = 1,2 - generation index of the lighter valence quark in 

the 5° meson, i.e. i = 2 chooses 5° and i = 1 chooses B^. 

Output delta_mb = ArriB^ for i = 1 

delta_mb = AniB^ for i = 2 

QCD related factors: 

coiiimon/meson_data/dmk,amk,epsk,fk,dmd,aind,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

Measured AM£^ dinb(l) = 3.01 ■ lO'^^ 

Measured AM^f dmb(2) = 1.2 ■ 10"" 

Measured width r^J gain_b(l) = 1.53 • lO-^^ 

Measured width T^^ gam-b(l) = 1-466 ■ 10"^^ 

coinmon/bx_4q/bk(5) ,bd(5),bb(2,5) , ainu_k , amu_d , amu_b 

BY^^ifiB) bb(l, 1) = bb(2, 1) = 0.87 

Bf^^{fiB) bb(l, 2) = bb(2, 2) = 0.8 

5fLL(/iB) bb(l, 3) = bb(2, 3) = 0.71 

Bl^iHB) bb(l, 4) = bb(2, 4) = 1.71 

B\^{^b) bb(l, 5) = bb(2, 5) = 1.16 

Renormalization scale ^b amu_b = 4.6 

common/ sm_4q/et a_c c,eta_ct,et a_t t , et a_b , bk_sm , bd_sm , bb_sm ( 2 ) 

BJ\Ib, bb_sm(l) = 1.22 

bb_sm(2) = 1.22 

r/fe eta_b = 0.55 

Details of calculations: Ref. [10] 



5.8 5^ — > Xg^ decay rate 

Both the SUSY contributions and the QCD corrections to the calculation of the B^ — )■ 
decay rate are quite complex. Their implementation in SUSY_FLAVOR is based on the SUSY 
loop calculations performed by the authors (not published in a general form) and on the 
QCD evolution published in [57] . There are no user-accessible QCD factors apart from the 
arguments of the bxg_nl routine. 

Routine double precision function bxg_nl(del , amiu_b) 

Input del - relative photon energy infrared cutoff scale, 

E^>{1- del)^^'^^-, < del< 1 

amiu_b - renormalization scale 
Output Br{B Xs-f). 

Details of calculations: General SUSY diagrams unpublished, QCD correc- 

tions based on [57] 
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6 Summary and Outlook 



We have presented SUSY_FLAVOR, a tool for calculating the set of important FCNC and CPV 
observables in the general i?-parity conserving MSSM. All implemented physical quantities 
(listed in Table [1]) can be calculated simultaneously for a given set of MSSM parameters. 
The calculations of the SUSY particle spectrum and flavor mixing matrices are performed 
exactly, so the code can be used for a completely general pattern of soft SUSY flavor 
violating terms and complex phases, without restrictions on the size of sfermion mass 
insertions. 

Besides complete routines for calculating the physical observables, SUSY_FLAVOR also pro- 
vides an extensive library of parton-level Green's functions and Wilson coefficients of many 
effective quark and lepton operators (see Table [3]). This set actually contains many more 
functions than are necessary to compute the quantities listed in Table [U These intermedi- 
ate building blocks can be used by SUSY_FLAVOR users to construct amplitudes for processes 
beyond those already fully implemented by dressing appropriate combinations of available 
form factors in QCD corrections and hadronic matrix elements, without repeating tedious 
one-loop SUSY calculations from scratch. For instance, the form factors implemented in 
SUSY_FLAVOR for the analysis of -B — >■ and Bd(s) — ^ decays [3l[l2] are sufficient to 
also calculate the B — )■ Kl'^l~ decay rate. 

SUSY_FLAVOR internally uses the conventions of Ref. [H], however in order to facilitate com- 
parison with other programs that analyze various sectors of MSSM, we have implemented 
an option to input parameters in the SLHA2 format [12]. 

SUSY_FLAVOR has been written in FORTRAN 77 and runs fairly quickly; it is capable of 
producing a reasonably wide-range scan over the MSSM parameters within hours or days 
on a typical personal computer. 

The SUSY_FLAVOR library is an open project. We want to gradually add more features in 
its future versions. In particular, we plan to: 

• add more observables in the 5-meson system, like the CP asymmetries in BB meson 
mixing and in B X^'y decay, as well as observables associated with B — Kl^l~ 
decay. 

• add observables for lepton flavor-violating processes like — )■ £^7, l'^ — )■ i^i^i'^ , and 
for the lepton anomalous magnetic moments, {g — 2)j 

• include quantities related to FCNCs in the top sector, like t — )■ cX with X = 
7, Z, g, H, in order to probe the flavor violation in up-squark mass matrices that 
are (almost) unconstrained to this moment. 

• implement full resummation of leading large tan /3 effects beyond the MFV scenario. 
With the increasing accuracy of experimental data on flavor and CP violation in rare 
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processes, it may eventually become possible to not only constrain the MSSM parameters, 
but also, if significant deviations from the SM predictions are found, to recover their actual 
values. For that multi-process analysis, such as the one performed by SUSY_FLAVOR, will be 
necessary. Therefore, we hope that SUSY_FLAVOR becomes an important tool that is useful 
not only to theorists working on MSSM but also to experimentalists fitting the MSSM onto 
forthcoming data from the Tevatron, LHC, and S-factories. 
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A Installation of the program 

The installation and execution of SUSY_FLAVOR is very simple. On Unix or Linux systems, 
just follow these steps : 



1. Download the code from http://www.fuw.edu.pl/susy_flavor and unpack it. 



2. Change directory into susy_f lavor. 

3. Edit Makefile and change F77 = gf ortran and FOPT = -0 -f no-automatic -Wall 
into your compiler name and options, respectively. 

4. Exit Makefile and type make (or gmake). 

5. If everything go through the code will ask you whether to read the input file susy_f lavor . in 
or to use the parameters defined inside the driver file. 

6. To run the code from now on just type ./sf lav. 

The authors tested SUSY_FLAVOR on Linux machines. With few straightforward modifi- 
cations the procedure describe above can be adapted to install program on other sys- 
tems. 
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The output of the program is displayed on the screen. In addition a file named mssm_data . txt 
is created. It contains information about the MSSM Lagrangian parameters and the tree- 
level mass spectrum corresponding to the input parameter set. A sample set of input param- 
eters and corresponding SUSY_FLAVOR output are listed in the following appendices. 

B Example of the SUSY_FLAVOR initialization sequence 

Below we present the contents of susy_f lavor . f , the master driver file for the SUSY_FLAVOR 
library. The driver program illustrates the correct initialization sequence for all relevant 
MSSM parameters (see SectionH]) and shows how to perform calls to the routines calculating 
physical observables (Section [5]). 

The driver file asks if the input parameters should be given directly inside the program or 
read from the default input file named susy_f lavor . in (in this case skipping the values 
given in the program). Defining the input parameters in the separate file is probably more 
straightforward, but the ability to initialize parameters from within the program could be 
more useful for performing mult i- dimensional scans over the MSSM parameter space. 

program susyjlavor 

implicit double precision (a-h,o-z) 

dimension sll(3) ,slr(3) ,sql(3) ,squ(3) ,sqd(3) 

double complex asl(3) ,asu(3) ,asd(3) 

double complex slmi_l(3) ,slmi_r(3) ,slmi_lr(3,3) 

double complex sqmi_l(3) , sdmi_r(3) , sumi_r(3) 

double complex sdmi_lr(3,3) , sumi_lr (3,3) 

double complex amg , amgg , amue 

common/ sf _cont/eps , indx (3,3) , iconv 

c decide if input parameters are read from file susy_flavor. in or defined inside the program 

write (*,' (a, $)') 'Read input from file susy_f lavor . in (no=l,yes=2)? ' 
read(*,*) input_type 
if (input_type . eq. 2) then 



call sflav_input 
goto 100 
end if 



Parameters read from file susy_fiavor. in 



c 



c 



Parameters defined inside the code. Start from input convention choice 

iconv = 1 ! SLHA2 input conventions 

iconv = 2 Jiep-ph/9511250^ input conventions 



c 



SM basic input initialization 



wmO = 80.398d0 



zmO = 91.1876d0 



M_Z 
! M_W 
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alphajz = 1/127. 934d0 / alpha_em(M_Z) 

call vpar_update(zmO,wmO,alpha_z) 



c QCD parameters 

alpha_s = 0.1172d0 
call lam_f it (alpha_s) 
call lam_f it jilo (alpha_s) 



/ alpha_s(M_Z) 

! fits Lambda-QCD at 3 loop level 
! fits Lambda_QCD at NLO level 



c CKM matrix initialization 

alam = 0.2258d0 / lambda 

apar = 0.808d0 / A 

rhobar = 0.177d0 / rho bar 

etabar = 0.360d0 / eta bar 

call ckm_wolf (alam, apar, rhobar, etabar) 

c Fermion mass initialization, input: MSbar running quark masses 

top_scale = 163. 2d0 

top = 163. 2d0 m_t(top_scale) 
bot_scale = 4.17d0 

bot = 4.17d0 m_b(bot_scale) 

call init_f erniion_sector (top,top_scale ,bot ,bot_scale) 

c Higgs sector parameters 

pm = 200 M_A 

tanbe =10 tan(beta) 

amue = (200.d0, lOO.dO) ! mu parameter 

call initJiiggs_sector (pm , tanbe , amue , ierr ) 

if (ierr.ne.O) stop 'negative tree level Higgs mass^?' 

c Gaugino sector parameters: if M1=0 set here then program uses 

amgg = (0.dO,0.dO) Ml (bino mass, complex) 

amg = (200.d0,0.d0) M2 (wino mass, complex) 

amglu = 3*abs(amg) M3 (gluino mass) 

call init_ino_sect or (amgg , amg , amglu , amue , tanbe , ierr) 
if (ierr.ne.O) write (*,*) '-ino mass below M_Z/2?' 



Slepton diagonal soft breaking parameters 

sll(l) = 300. dO left selectron mass scale 



sll(2) = 300. dO 

sll(3) = 300. dO 

slr(l) = 300. dO 

sir (2) = 300. dO 

sir (3) = 300. dO 



left smuon mass scale 
left stau mass scale 
right selectron mass scale 
right smuon mass scale 
right stau mass scale 



Dimensionless (normalized to masses) slepton diagonal LR mixing 
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asl(l) = (1 .dO.O.dO) / 1st generation 

asl(2) = (1 .dO,0.dO) ! 2nd generation 

asl(3) = (l.dO.O.dO) ! 3rd generation 

Slepton LL and RR mass insertions (hermitian matrices, only upper part given) 
slmi_x(l),slmi_x(2), slmi_x(3) are 12,23,31 entry, respectively 
do i=l,3 

slmi_l(i) = (0.dO,0.dO) ! slepton LL mass insertion 

slmi_r(i) = (0.dO,0.dO) / slepton RR mass insertion 

end do 

sliiii_l(2) = (2 . d-2 , 1 . d-2) / example, non-vanishing LL 23 entry 

Slepton LR mass insertions, non-hermitian in general 
do i=l,3 
do j=l,3 

slini_lr(i, j) = (0.dO,0.dO) / slepton LR ij mass insertion 
end do 
end do 

Calculate slepton physical masses and mixing angles 

call init_slepton_sector (sll , sir , asl , ierr , slmi J. , slmi_r , slmi_lr) 
if (ierr.ne.O) stop 'negative tree level slepton mass^?' 



Squark diagonal soft breaking parameters 



sql(l) 


= 500 


dO 


/ left squark mass, 1st generation 


sql(2) 


= 500 


dO 


! left squark mass, 2nd generation 


sql(3) 


= 400 


dO 


! left squark mass, 3rd generation 


sqd(l) 


= 550 


dO 


! right down squark mass 


sqd(2) 


= 550 


dO 


! right strange squark mass 


sqd(3) 


= 300 


dO 


! right sbottom mass 


squ(l) 


= 450 


dO 


! right up squark mass 


squ(2) 


= 450 


dO 


! right charm squark mass 


squ(3) 


= 200 


dO 


! right stop mass 



Dimensionless (normalized to masses 
asd(l) = (l.dO.O.dO) 
(l.dO.O.dO) 
(l.dO.O.dO) 
(l.dO.O.dO) 
(l.dO.O.dO) 
(l.dO.O.dO) 



asd(2) 
asd(3) 
asu(l) 
asu(2) 
asu(3) 



') squark diagonal LR mixing 
! down squark LR mixing, 1st generation 
! down squark LR mixing, 2nd generation 
! down squark LR mixing, 3rd generation 
up squark LR mixing, 1st generation 
up squark LR mixing, 2nd generation 
up squark LR mixing, 3rd generation 
Squark LL and RR mass insertions (hermitian matrices, only upper part given) 
sqmiJ(l),sqmiJ(2), sqmiJ(3) are 12,23,31 entry, respectively, etc. 
do i=l,3 

sqmi_l(i) = (O.dO.O.dO) ! squark LL mass insertion 

sunii_r(i) = (O.dO.O.dO) ! up- squark RR mass insertion 

sdmi_r(i) = (O.dO.O.dO) down-squark RR mass insertion 

end do 
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sqmi_l(2) = (2 . d-2 , -1 . d-2) / example, non-vanishing LL 23 entry 

c Squark LR mass insertions, non-hermitian in general 

do i=l,3 
do j=l,3 

suini_lr(i, j) = (0.dO,0.dO) ! up-squark LR ij mass insertion 
sdmi_lr(i, j) = (0.dO,0.dO) / down-squark LR ij mass insertion 
end do 
end do 

c Calculate squark physical masses and mixing angles 

call init_squark_sector ( sql , squ , sqd , asu , asd , ierr , sqmi J. , sumi_r , 
$ sdini_r,sumi_lr,sdmi_lr) 
if (ierr.ne.O) stop 'negative tree level squark mass^?' 

c reset status of physical Higgs mass after parameter changes 

call reset_phys_data 
c Neutral CP-even Higgs masses in the 1-loop Effective Potential Approximation, 

c Only real mu, A_t, AJ) allowed - replace x-^ahs(x) 

call f corr_EPA(tanbe,pm,top,abs(aniue) ,sql(3) ,sqd(3) ,squ(3) , 
$ abs (asd(3) ) , abs (asu(3) ) , ierr) 

if (ierr.ne.O) stop 'negative 1-loop EPA CP-even Higgs mass^?' 

100 continue End of input section !!! 

c Control output: Lagrangian parameters and tree level masses written on file mssm_data.txt 

ifl = 1 output file number 

open ( if 1 , f ile= ' mssm_data . txt ' , status= ' unknown ' ) 
call print_MSSM_par (if 1) / Lagrangian parameters 

call print_MSSM_niasses (if 1) / tree level physical masses 
close (if 1) 

c Results for implemented rare decays: 

write(*, 99) 'Electric dipole moments:' 
write(*, 99) 'Electron EDM = ',edm_l(l) 
write(*,99) 'Muon EDM = ',edm_l(2) 
write (*, 99) 'Tau EDM = ',edm_l(3) 
write (*, 99) 'Neutron EDM = ',edm_n() 

write (*, 99) 'Neutrino K decays:' 

call k_pivv(br_kO,br_kp) 

write(*,99) 'BR(k5, pi° vv) = ' ,br_kO 

write(*,99) 'BR(K+ pi+ vv) = ' ,br_kp 

write (*, 99) 'Leptonic B decays:' 

write (*, 99) 'BR(B_d mu+ mu") = ' ,b_ll(3, 1 ,2,2) 
write(*,99) 'BR(B^ mu+ mu~) = ' ,b_ll(3,2,2,2) 
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write(*,99) 'B— 7- X_s photon decay:' 

c Physical quantities for BR(B^X_s g) calculation 

delb = 0.99d0 / Photon energy infrared cutoff 

amiu_b= 4.8d0 / Renormalization scale miu-b 

write(*,99) 'BR(B X_S gamma) = ' ,bxg_nl(delb,amiu_b) 

write (*, 99) 'KK mixing:' 
call dd_kaon(epsJs,delta_mk) 
write (*, 99) 'eps_K = ',epsjc 
write (*, 99) 'Delta m_K = ' ,deltajnk 

write(*,99) 'DD mixing:' 
call uu_dmeson(delta_md) 
write (*, 99) 'Delta m_D = ' .deltajnd 

write(*,99) 'BB mixing:' 

call dd_bmeson(l ,deltajnbd) 

write(*, 99) 'Delta m_B_d = ' .deltajnbd 

call dd_bmeson(2,deltajnbs) 

write (*, 99) 'Delta m_B_s = ' ,delta_nibs 

99 f ormat(a, Ipell .4) 
end 

C Example of SUSY_FLAVOR input file 

By default, the driver file susy_f lavor . f reads input parameters from the file susy_f lavor . in. 
Below we provide an example input file defining a set of parameters equivalent to those in 
the driver file presented in Appendix [Bl 

# Example input of SUSY_FLAVOR in Les Houches Accord-like format 
# 

# CAUTION: users can modify numerical data in this file but they 

# should not remove existing data lines within blocks SMINPUTS, 

# VCKMIN, EXTPAR, MSL2IN, MSE2IN, MSQ2IN, MSU2IN, MSD2IN, TEIN, TUIN, 

# TDIN, IMMSL2IN, IMMSE2IN, IMMSQ2IN, IMMSU2IN, IMMSD2IN, IMTEIN, 

# IMTUIN, IMTDIN. New data lines in each block can be added but only 

# after the already defined ones. Also, comment-only lines starting 

# from # as a first character can be added only just after or before 

# Block XXX statements, i.e. not between data lines. Otherwise 

# SUSY_FLAVOR input routine sflav_input will denounce input file as 
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# corrupted or read incorrect values . 
# 

# Full new data blocks can be added, sflav_input will ignore them. 
# 

Block MODSEL # Select model 

10 # General MSSM 

3 # MSSM particle content 

4 # R-parity conserving MSSM 

5 2 # CP violated 

6 3 # Lepton and quark flavor violated 
Block SOFTINP # Choose convention for the soft terms 

# convention = 1: 

# sfermion input parameters in SLHA2 conventions 

# convention = 2: 

# sfermion input parameters in conventions of hep-ph/9511250 

# input_type = 1 : 

# sfermion off-diagonal terms given as dimensionless mass insertions 

# LR diagonal terms given as dimensionless parameters 

# input_type = 2 : 

# sfermion soft terms given as absolute values 

# See comment in Blocks MSXIN2, TXIN below 
1 



2 

2 1 
Block SMINPUTS 



# sfermion convention, SLHA2 or hep-ph/9511250 



# input_type (dimension of soft mass entries) 

# Standard Model inputs 



1 


1 


279340000e+02 


# 


alpha" 1 SM MSbar(MZ) 


3 


1 


172000000e-01 


# 


alpha_s(MZ) SM MSbar 


4 


9 


118760000e+01 


# 


MZ(pole) 


5 


4 


170000000e+00 


# 


mb(mb) SM MSbar 


6 


1 


632000000e+02 


# 


mtop(mt) SM MSbar 


7 


1 


777000000e+00 


# 


mtau(pole) 


11 


5 


llOOOOOOOe-04 


# 


me (pole) 


13 


1 


056590000e-01 


# 


mmu(pole) 


21 


7 


OOOOOOOOOe-03 


# 


md(2 GeV) MSbar 


22 


4 


OOOOOOOOOe-03 


# 


mu(2 GeV) MSbar 


23 


1 


lOOOOOOOOe-01 


# 


ms(2 GeV) MSbar 


24 


1 


279000000e+00 


# 


mc(mc) MSbar 


30 


8 


039800000e+01 


# 


MW (pole) , not standard 


:k VCKMIN 


# 


CKM matrix 


1 


2 


258000000e-01 


# 


lambda 


2 


8 


080000000e-01 


# 


A 


3 


1 


770000000e-01 


# 


rho bar 


4 


3 


600000000e-01 


# 


eta bar 



Block EXTPAR 

1 O.OOOOOOOOOe+02 



# non-minimal input parameters , real part 

# Re(ml), U(l) gaugino mass 
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2 2.000000000e+02 

3 6.000000000e+02 
23 2.000000000e+02 

25 l.OOOOOOOOOe+01 

26 2.0000000006+02 
Block IMEXTPAR 

1 O.OOOOOOOOOe+00 

2 O.OOOOOOOOOe+00 
23 l.OOOOOOOOOe+02 



# Re(m2), SU(2) gaugino mass 

# m3, SU(3) gaugino mass 

# Re(mu) 

# tan(beta) 

# MA 

# non-minimal input parameters, imaginary part 

# Im(ml) , U(l) gaugino mass 

# Im(m2) , SU(2) gaugino mass 

# Im(mu) 



# if abs(ml) = SUSY_FLAVOR uses ml=5/3 s_W/c_W m2 
# 

# Soft sfermion mass matrices 
# 

# Of f -diagonal entries may be given as absolute entries or as 

# dimensionless mass insertions - then real off-diagonal entries of 

# SLHA2 blocks are calculated by SUSY_FLAVDUR as 

# M^CI.J) = (mass insertion) (I , J) sqrtCM^CD M^CJ.J)) 

# (see comments at the top of subroutine sflav_input) 
# 

# Below we give an example of dimensionless off-diagonal entries 
# 



Block 


MSL2IN 


# 


left soft slepton mass matrix, real part 


1 


1 9.0000000006+04 


# 


Left slepton diagonal mass^, 1st generation 


2 


2 9.0000000006+04 


# 


Left slepton diagonal mass^, 2nd generation 


3 


3 9.0000000006+04 


# 


Left slepton diagonal mass^, 3rd generation 


1 


2 O.OOOOOOOOOe+OO 


# 


Dimensionless left slepton mass insertion 12 


2 


3 2.0000000006-02 


# 


Dimensionless left slepton mass insertion 23 


1 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless left slepton mass insertion 13 


Block 


IMMSL2IN 


# 


left soft slepton mass matrix, imaginary part 


1 


2 O.OOOOOOOOOe+OO 


# 


Dimensionless left slepton mass insertion 12 


2 


3 1.0000000006-02 


# 


Dimensionless left slepton mass insertion 23 


1 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless left slepton mass insertion 13 


Block 


MSE2IN 


# 


right soft slepton mass matrix, real part 


1 


1 9.0000000006+04 


# 


Right selectron diagonal mass^ 


2 


2 9.0000000006+04 


# 


Right smuon diagonal mass^ 


3 


3 9.0000000006+04 


# 


Right stau diagonal mass^ 


1 


2 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 12 


2 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 23 


1 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 13 


Block 


IMMSE2IN 


# 


right soft slepton mass matrix, imaginary part 


1 


2 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 12 


2 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 23 


1 


3 O.OOOOOOOOOe+OO 


# 


Dimensionless right slepton mass insertion 13 


Block 


MSQ2IN 


# 


left soft squark mass matrix, real part 
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1 1 2.500000000e+05 

2 2 2.500000000e+05 

3 3 1.600000000e+05 

1 2 O.OOOOOOOOOe+00 

2 3 2.000000000e-02 
1 3 O.OOOOOOOOOe+00 

Block IMMSQ2IN 

1 2 O.OOOOOOOOOe+00 

2 3 -l.OOOOOOOOOe-02 
1 3 O.OOOOOOOOOe+00 

Block MSU2IN 

1 1 2.025000000e+05 

2 2 2.025000000e+05 

3 3 4.000000000e+04 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 

Block IMMSU2IN 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 



# Left squark diagonal mass , 1st generation 

# Left squark diagonal mass^, 2nd generation 

# Left squark diagonal mass^, 3rd generation 

# Dimensionless left squark mass insertion 12 

# Dimensionless left squark mass insertion 23 

# Dimensionless left squark mass insertion 13 

# left soft squark mass matrix, imaginary part 

# Dimensionless left squark mass insertion 12 

# Dimensionless left squark mass insertion 23 

# Dimensionless left squark mass insertion 13 

# right soft up-squark mass matrix, real part 

# Right u-squark diagonal mass^ 

# Right c-squark diagonal mass^ 

# Right stop diagonal mass^ 

# Dimensionless right up-squark mass insertion 12 

# Dimensionless right up-squark mass insertion 23 

# Dimensionless right up-squark mass insertion 13 

# right soft up-squark mass matrix, imaginary part 

# Dimensionless right up-squark mass insertion 12 

# Dimensionless right up-squark mass insertion 23 

# Dimensionless right up-squark mass insertion 13 



Block MSD2IN # right soft down-squark mass matrix, real part 



1 1 3.025000000e+05 

2 2 3.025000000e+05 

3 3 9.000000000e+04 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 

Block IMMSD2IN 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 



# Right d-squark diagonal mass 

# Right s-squark diagonal mass^ 

# Right sbottom diagonal mass^ 

# Dimensionless right down-squark mass insertion 12 

# Dimensionless right down-squark mass insertion 23 

# Dimensionless right down-squark mass insertion 13 

# right soft down-squark mass matrix, imaginary part 

# Dimensionless right down-squark mass insertion 12 

# Dimensionless right down-squark mass insertion 23 

# Dimensionless right down-squark mass insertion 13 



# 

# Soft sfermion trilinear mixing matrices 
# 

# LR mixing parameters can be given as absolute entries or as 

# dimensionless diagonal A-terms and dimensionless ff-diagonal mass 

# insertions - see comments at the top of subroutine sflav_input 
# 

# Below we give an example of dimensionless A terms. 
# 

# Diagonal entries below are dimensionless "A parameters" 

# Diagonal entries of SLHA2 LR blocks are calculated by SUSY_FLAVOUR as 

# TL(I,I) = AL(I,I) Yukawa_L(I) sqrt(ML2(I,I)*ME2(I,I)) 
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# TU(I,I) = AU(I,I) Yukawa_U(I) sqrt (MQ^CI , I) *MU2 (I , I) ) 

# TD(I,I) = AD(I,I) Yukawa_D(I) sqrt (MQ^CI , I) *MD2 (I , I) ) 
# 

# Of f -diagonal entries are dimensionless "mass insertions" 

# Off-diagonal entries of SLHA2 LR blocks are calculated by SUSY_FLAVOUR as 
# 

# TL(I,J) = AL(I,J) sqrt(2 ML^ (I , I) *ME2 ( J , J) ) /vl 

# TU(I,J) = AU(I,J) sqrt(2 MQ^ (I , I) *MU2 ( J , J) ) /v2 

# TD(I,J) = AD(I,J) sqrt(2 MQ^ (I , I) *MD2 ( J , J) ) /vl 
# 

Block TEIN # slepton trilinear mixing, dimensionless, real part 



1 


1 


1. 


,000000000e+00 


# 


Diagonal AL term, 1st generation 




2 


2 


1. 


,000000000e+00 


# 


Diagonal AL term, 2nd generation 




3 


3 


1. 


,000000000e+00 


# 


Diagonal AL term, 3rd generation 




1 


2 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 12 




2 


1 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 21 




2 


3 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 23 




3 


2 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 32 




1 


3 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 13 




3 


1 


0. 


OOOOOOOOOe+OO 


# 


Slepton LR mass insertion 31 




ck 


IMTEIN 


# 


slepton trilinear mixing, dimensionless. 


imag. 


1 


1 


0. 


,000000000e+00 


# 


Diagonal AL term, 1st generation 




2 


2 


0. 


,000000000e+00 


# 


Diagonal AL term, 2nd generation 




3 


3 


0. 


OOOOOOOOOe+OO 


# 


Diagonal AL term, 3rd generation 




1 


2 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 12 




2 


1 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 21 




2 


3 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 23 




3 


2 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 32 




1 


3 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 13 




3 


1 


0. 


,000000000e+00 


# 


Slepton LR mass insertion 31 




ck 


TUIN 


# 


up-squark trilinear mixing, dimensionless 


, real 


1 


1 


1. 


,000000000e+00 


# 


Diagonal AU term, 1st generation 




2 


2 


1. 


,000000000e+00 


# 


Diagonal AU term, 2nd generation 




3 


3 


1. 


,000000000e+00 


# 


Diagonal AU term, 3rd generation 




1 


2 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 12 




2 


1 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 21 




2 


3 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 23 




3 


2 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 32 




1 


3 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 13 




3 


1 


0. 


,000000000e+00 


# 


Up-squark LR mass insertion 31 




ck 


IMTUIN 


# 


up-squark trilinear mixing, dimensionless 


, imag 


1 


1 


0. 


,000000000e+00 


# 


Diagonal AU term, 1st generation 




2 


2 


0. 


OOOOOOOOOe+OO 


# 


Diagonal AU term, 2nd generation 




3 


3 


0. 


,000000000e+00 


# 


Diagonal AU term, 3rd generation 





part 



part 
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1 


2 


0. 


OOOOOOOOOe+00 


# 


2 


1 


0. 


OOOOOOOOOe+00 


# 


2 


3 


0. 


OOOOOOOOOe+00 


# 


3 


2 


0. 


OOOOOOOOOe+00 


# 


1 


3 


0. 


OOOOOOOOOe+00 


# 


3 


1 


0. 


OOOOOOOOOe+00 


# 


ck 


TDIN 


# 


1 


1 


1. 


OOOOOOOOOe+00 


# 


2 


2 


1. 


OOOOOOOOOe+00 


# 


3 


3 


1. 


OOOOOOOOOe+00 


# 


1 


2 


0. 


OOOOOOOOOe+00 


# 


2 


1 


0. 


OOOOOOOOOe+OO 


# 


2 


3 


0. 


OOOOOOOOOe+00 


# 


3 


2 


0. 


OOOOOOOOOe+00 


# 


1 


3 


0. 


OOOOOOOOOe+00 


# 


3 


1 


0. 


OOOOOOOOOe+OO 


# 


ck 


IMTDIN 


# 


1 


1 


0. 


, OOOOOOOOOe+OO 


# 


2 


2 


0. 


, OOOOOOOOOe+OO 


# 


3 


3 


0. 


, OOOOOOOOOe+OO 


# 


1 


2 


0. 


OOOOOOOOOe+OO 


# 


2 


1 


0. 


OOOOOOOOOe+OO 


# 


2 


3 


0. 


OOOOOOOOOe+OO 


# 


3 


2 


0. 


OOOOOOOOOe+OO 


# 


1 


3 


0. 


OOOOOOOOOe+OO 


# 


3 


1 


0. 


OOOOOOOOOe+OO 


# 



D Example of SUSY_FLAVOR output 

The parameters defined inside the driver program in Appendix [B] and in the input file hsted 
in Appendix [Q should produce identical output. We enclose it here so that SUSY_FLAVOR 
users can check that the program gives the same result on their own computers and FOR- 
TRAN compilers. 

The driver file susy_f lavor . f writes the MSSM Lagrangian parameters and tree-level 
particle masses to the file mssni_data.txt. For the parameters defined in Appendices iBl 
and O one has: 

******* MSSM Lagrangian parameters ******* 
QED coupling l/alpha_em(M_Z) = 1.2793E+02 

Weinberg angle s_W^ = 2.2265E-01 

Z boson mass = 9.1188E+01 
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W boson mass 

QCD coupling alpha_s (M_Z) 

Higgs mixing parameter mu (complex) 
Higgs soft mixing parameter m^2 
Higgs soft masses >^%2 

U(l) gaugino mass (complex) 
SU(2) gaugino mass (complex) 
SU(3) gaugino mass (real) 

Left slepton mass matrix, real part 
9.00009E+04 O.OOOOOE+00 . OOOOOE+00 

O.OOOOOE+00 9.00018E+04 1 . 80004E+03 

O.OOOOOE+00 1.80004E+03 9.00027E+04 

Left slepton mass matrix, imaginary part: 
O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

-O.OOOOOE+00 O.OOOOOE+00 9 . 00022E+02 

-O.OOOOOE+00 -9.00022E+02 O.OOOOOE+00 

Right slepton mass matrix, real part: 

8.99991E+04 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 8.99982E+04 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 8 . 99973E+04 

Right slepton mass matrix, imaginary part: 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

Slepton LR mixing matrix, real part: 

-9.00010E-03 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 -1.86094E+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 -3.12978E+01 

Slepton LR mixing matrix, imaginary part: 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

Left squark mass matrix, real part: 
2.50003E+05 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 2.50005E+05 4.00010E+03 

O.OOOOOE+00 4.00010E+03 1 . 60005E+05 

Left squark mass matrix, imaginary part: 
O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 



= 8.0398E+01 
= 1.1720E-01 

= 2.0000E+02 
= -3.9604E+03 
= -6.3208E+03 

= 9.5472E+01 
= 2.0000E+02 
= 6.0000E+02 



l.OOOOE+02 

-5.3679E+04 

O.OOOOE+00 
O.OOOOE+00 
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O.OOOOOE+00 O.OOOOOE+00 -2 . 00005E+03 

O.OOOOOE+00 2.00005E+03 O.OOOOOE+00 



Right up-squark 

2.02498E+05 

O.OOOOOE+00 

O.OOOOOE+00 

Right up-squark 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 



mass matrix, 

O.OOOOOE+00 

2.02496E+05 

O.OOOOOE+00 

mass matrix, 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 



real part : 

O.OOOOOE+00 
O.OOOOOE+00 
3.99988E+04 

imaginary part : 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 



Right down-squark mass matrix, real part: 

3.02497E+05 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 3.02494E+05 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 8 . 99973E+04 

Right down-squark mass matrix, imaginary part: 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 



Up-squark LR 

6.18168E-03 

O.OOOOOE+00 

O.OOOOOE+00 

Up-squark LR 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 



mixing matrix, 
O.OOOOOE+00 
1.79332E+00 
O.OOOOOE+00 

mixing matrix, 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 



real part : 

O.OOOOOE+00 
O.OOOOOE+00 
2.71002E+02 

imaginary part : 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 



Down-squark LR 

-1.19597E-01 

O.OOOOOE+00 

O.OOOOOE+00 

Down-squark LR 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 



mixing matrix, 
O.OOOOOE+00 
-1.87938E+00 
O.OOOOOE+00 

mixing matrix, 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 



real part : 
O.OOOOOE+00 
O.OOOOOE+00 
-5.56628E+01 

imaginary part : 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 



****** Particle masses 
* Fermion masses ** 
Charged lepton masses 
Running u quark masses 
Running d quark masses 



in GeV : ******* 

5.110E-04 
at nit scale 2.220E-03 
at nit scale 3.885E-03 



1.057E-01 1.777E+00 
6.440E-01 1.632E+02 
6.104E-02 2.737E+00 
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* Higgs masses ** 

Tree level (H,h,A,H+) : 

1-loop, EPA approximation (H,h) 



2.010E+02 8.893E+01 2.000E+02 2.156E+02 
2.005E+02 1.137E+02 



* Tree level SUSY masses ** 



Sneutrino masses 


2. 


, 897E+02 


2, 


.931E+02 


2. 


,965E+02 












Slepton masses 


2. 


,953E+02 


3, 


.028E+02 


3. 


, 030E+02 


3, 


.037E+02 


3. 


.038E+02 


3 


U squark masses 


2. 


. 178E+02 


4, 


.486E+02 


4. 


,487E+02 


4, 


,489E+02 


4. 


.971E+02 


4 


D squark masses 


2. 


.999E+02 


4, 


.049E+02 


5. 


, 035E+02 


5, 


.037E+02 


5. 


.505E+02 


5 


Chargino masses 


1, 


,552E+02 


2, 


.808E+02 
















Neutralino masses 


8. 


.865E+01 


1, 


.584E+02 


2. 


,322E+02 


2, 


.808E+02 








Gluino mass 


6. 


.OOOE+02 





















The output for the physical observables is printed on the standard output, usually the 
computer screen of the console. It should look like: 



Electric dipole moments: 
Electron EDM 
Muon EDM 
Tau EDM 
Neutron EDM 



= 4.7256E-25 
= 9.7726E-23 
= 1.6425E-21 
= 5.9331E-24 



Neutrino K decays: 
BR(KX° -> pi° vv) 
BR(K+ -> pi+ vv) 

Leptonic B decays: 
BR(B_d -> mu+ mu~) 
BR(B_s -> mu+ mu~) 



= 2.8555E-11 
= 7.3932E-11 



= 1.2012E-10 
= 4.7395E-09 



B -> X_s photon decay: 
BR(B -> X_s gamma) 



= 2.5756E-04 



KK mixing: 
eps_K 
Delta m_K 



= 2.3366E-03 
= 2.4362E-15 



DD mixing: 
Delta m_D 



= 1.6656E-17 



BB mixing: 
Delta m_B_d 
Delta m_B_s 



= 3.6999E-13 
= 1.3242E-11 
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PROGRAM SUMMARY 

Manuscript Title: SUSY_FLAVOR: a computational tool for FCNC and CP-violating processes in 
the MSSM 

Authors: J. Rosiek, P. H. Chankowski, A. Dedes, S. Jager, P. Tanedo 
Program Title: SUSY_FLAVDR 

Journal Reference: 
Catalogue identifier: 
Licensing provisions: None 
Programming language: Fortran 77 
Operating system: Any, tested on Linux 

Keywords: Supersymmetry, K physics, B physics, rare decays, CP-violation 
PACS: 12.60. Jv, 13.20.He 

Classification: 11.6 Phenomenological and Empirical Models and Theories 
Nature of problem: 

Predicting CP-violating observables, meson mixing parameters and branching ratios for set of rare 
processes in the general R-parity conserving MSSM. 
Solution method: 

We use standard quantum theoretical methods to calculate Wilson coefficients in MSSM and at 
one loop including QCD corrections at higher orders when this is necessary and possible. The 
input parameters can be read from an external file in SLHA format. 
Restrictions: 

The results apply only to the case of MSSM with R-parity conservation. 
Unusual features: 
Running time: 

For single parameter set approximately Is in double precision on a PowerBook Mac G4 
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